Sub-picosecond compression by velocity bunching in a photo-injector 
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Abstract 

We present an experimental evidence of a bunch compression scheme that uses a traveling wave accelerating 
structure as a compressor. The bunch length issued from a laser-driven radio-frequency electron source was 
compressed by a factor >3 using an S-band traveling wave structure located immediately downstream from 
the electron source. Experimental data are found to be in good agreement with particle tracking simulations. 
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I. INTRODUCTION 



In the recent years there has been an increasing demand on ultrashort electron bunches to drive 
short-wavelength free-electron lasers and study novel accelerating techniques such as plasma-based 
accelerators [jl], g]. Short bunches are commonly obtained by magnetic compression. In this latter 
scheme, the bunch is compressed using a series of dipoles arranged in a chicane configuration such 
to introduce an energy-dependent pathlength. Therefore an electron bunch having the proper 
time-energy correlation can be shortened in the chicane. However, problems inherent to magnetic 
compression such as momentum spread and transverse emittance dilution due to the bunch self- 
interaction via coherent synchrotron radiation [||| has brought back the idea of bunching the beam 
with radio- frequency (rf) structures ||. 

It was recently proposed to incorporate the latter method (henceforth named velocity bunching) 
into the next photo-injector designs ||. The velocity bunching relies on the phase slippage between 
the electrons and the rf-wave that occurs during the acceleration of non ultra-relativistic electrons. 
In this paper after presenting a brief analysis of the velocity bunching scheme, we report on its 
exploration at the deep ultraviolet free-electron laser (DUV-FEL) facility of Brookhaven National 
Laboratory (BNL). The measurements are compared with numerical simulations performed with 
the computer program ASTRA J/J. 

II. ANALYSIS OF THE VELOCITY BUNCHING TECHNIQUE 

In this section we elaborate a simple model that describes how the velocity bunching works. A 
more detailed discussion is given in Reference Q. 

An electron in an rf traveling wave accelerating structure experiences the longitudinal electric field: 



where E Q is the peak field, k the rf wavenumber and ip the injection phase of the electron with 
respect to the rf wave. Let ip(z, t) = cut — kz + tp be the relative phase of the electron w.r.t the 
wave. The evolution of ip(t, z) can be expressed as a function of z solely: 



E z (z, t) = E Q sui{ujt — kz + V> ) 



(1) 





Introducing the parameter a = fc eg ° a , we write for the energy gradient [|| : 

Y = ak sua (V0- (3) 

The system of coupled differential equations (H) and (|3|) with the initial conditions j(z = 0) = 7 Q 
and ip(z = Q) = ip describe the longitudinal motion of an electron in the rf structure. Such a 
system is solved using the variable separation technique to yield: 



a cos ip + j- \/7 2 - l = C. (4) 
Or, expliciting ^ as a function of 7: 



^(7) = arccos [ C 7jVg_ [] _ (5) 



o 



Here the constant of integration is set by the initial conditions of the problem|Lq]: C = acosip a + 
7o — ^7^ — 1. The latter equation gives insights on the underlying mechanism that provides 
compression. In order to get a simpler model, we consider the limit: ipoo = h m 7^oo^(7) = 
arccos ^cos(^ ) + 2^^) i we nave assumed j is larger than unit and did the approximation 7 Q — 
\/"io — 1 — l/(27 ). After differentiation of Eq. |^, given an initial phase dtp and energy d^ extents 
we have for the final phase extent: 

sin(V'o) , , , 1 , {R . 

dipoc = -T—-—-rdip + - — . d7 Q . (6) 

sm(V'oo) 2a7^sm(V'oo) 

Hence depending upon the incoming energy and phase extents, the phase of injection in the rf 
structure ip can be tuned to minimize the phase extent after extraction i.e. to ideally (under 
single-particle dynamics) make dipoo — > 0. We note that there are two contributions to cfa/'ocv the 
first term dip ^1 dip comes from the phase slippage (the injection and extraction phases are generally 
different). The second term dipoo/d'fo is the contribution coming from the initial energy spread. 
To illustrate the compression mechanism we consider a two macro-particles model. In Figure [l] we 
present results obtained by numerically integrating the equation of motion for two non-interacting 
macro-particles injected into a 3 m long traveling wave structure. Given the incoming phase Aip a 
and energy Ajo spreads between the two macro-particles, and the accelerating gradient of the 
structure (taken to be 20 MV/m), we can optimize the injection phase to minimize the bunch 
length at the structure exit. 



III. EXPERIMENTAL RESULTS 



The measurement was carried out at the DUV-FEL facility of Brookhaven national laboratory [11]. 
A block diagram of the linear accelerator is given in Fig. |2[ The electron bunches of ~4 MeV, 
generated by a laser-driven rf electron source, are accelerated by a series of four linac sections. The 
linac sections consist of 2.856 GHz traveling wave structures operating on the 2n/3 accelerating 
mode. The structures are approximately 3 m long and can operate with an average accelerating 
voltage up to 20 MV/m. Nominally the bunch is shortened using a magnetic bunch compressor 
chicane located between the second and third linac sections. In this latter case, the linac sections 
LI, L3, L4 are ran on-crest while the linac L2 is operated to impart the proper time-energy 
correlation along the bunch to enable compression as the beam pass through the magnetic chicane. 



To investigate the velocity bunching scheme, the linac section LI was used as a buncher: its 
phase was varied and, for each phase setting, the section L2 was properly phased to maximize the 
beam energy with sections L3 and L4 turned off. The magnetic bunch compressor was turned off 
during the measurement. The nominal settings for the different rf and photo-cathode drive-laser 
parameters are gathered in Table |. 

parameter value units 

laser injection phase 40 ± 5 rf-deg 

laser radius on cathode 0.75 ± 0.1 mm 
laser rms length 1.15 ± 0.1 ps 

E-peak on cathode 83 ± 1 MV/m 

LI average accelerating field 10.5 ± 0.1 MV/m 

L2 average accelerating field 13.2 ± 0.1 MV/m 

TABLE I: Nominal settings for the rf-gun, accelerating sections, and the photo-cathode drive-laser. The 
values have been directly measured or inferred from the beam properties. 

The measurements of bunch length that follow are compared with simulations performed with the 
program ASTRA 0. ASTRA is a macro-particle tracking code based on a rotational symmetric 
space charge algorithm. It incorporates a detailed model for the traveling wave accelerating 
structure g, [J. To perform the simulations we used the parameters values of Table |. The laser 



transverse distribution was modeled by a radially uniform transverse distribution with 0.75 mm 
radius, and the time profile, measured using a single shot cross-correlation technique, was directly 
loaded into the simulations. 



Both time- and frequency-domain techniques were used to characterize the bunching process as the 
phase of the linac LI was varied. 



The time-domain charge density was directly measured using the so-called zero-phasing method [12, 



13]. In the present case, we use the linac section L3 to cancel the incoming time-energy correlation, 
and operate the linac L4 at zero-crossing to introduce a linear time-dependent energy chirp along 
the bunch (we have investigated both zero-crossing points). The bunch is then directed to a beam 
viewer ("YaG monitor" in Fig. downstream from a 72° angle spectrometer. The viewer, located at 
a dispersion (horizontal) of ij = 907 mm, allows the measurement of the bunch energy distribution. 
Unlike in Reference |l2| , the longitudinal phase space of beams issued from an rf electron source 
is not perfectly linear: because of the longitudinal space charge forces, the phase space generally 
has a third order distortion fll0"|| . To analyze the impact of such a distortion on our bunch length 
measurement method, it is interesting to consider the Gaussian normalized longitudinal phase space 
(s, 8) density: 

^/ o 1 ( (<* " his - h 3 s 3 ) 2 \ ( s 2 \ 

Here <j s and as are the bunch rms length and rms uncorrelated fractional momentum spread 
and h\, /13 are constants that quantify the linear and third order correlations of the longitudinal 
phase space. The zero-phasing measurement can then be analyzed in term of a sequence of 
numerical calculation based on Eq. [7J by computing and comparing the time and fractional 
momentum spread projections associated to 7^(5,(5 + C Q x s). The constant C Q depends on the 
incoming beam energy E Q , the accelerating voltage of the zero-phased linac section, the rf wavenum- 
ber k r f, and dispersion Q|: C Q = ± ^ y E ° k - the ± sign reflects the two possible zero-crossing points. 



An example of such a calculation is presented in Fig. [| To generate the presented data 
we started with a longitudinal phase space which has a third order distortion but no linear 
correlation (as it should be downstream from linac L3) we then set the constant C a to have 
a full-width fractional momentum spread of approximately 1.5 % similar to the value imposed 



by the finite size of the viewer (diameter ~15 mm) used for the measurement of the bunch 
energy distribution. The Figure ^ demonstrates the impact of the third order distortion in 
the longitudinal phase space: depending on the chosen zero-crossing phase, it contributes to 
an elongation or a contraction of the measured time profile compared to the real profile. For 
the rms bunch length measurements reported hereafter we computed the average bunch length 
measured for the two zero-crossing points and reported the difference as an error bar. For the 
reported bunch profiles we use the bunch profile corresponding to the case when the phase space has 
no fold over. Hence we expect the bunch time-profile reported hereafter to be longer than in reality. 

As the phase of the linac section LI was varied and L2 tuned to maximize the energy gain, the beam 
energy was measured. The so-obtained energy variation versus the phase of the linac LI is compared 
with the simulations for the nominal operating point (see Table |) in Fig. |] and the corresponding 
plot for the bunch length is shown in Fig. || As predicted, we observed that operating the linac at 
lower phases (thereby giving the bunch head less energy than the tail) provides some compression. 
The parametric dependence of the rms bunch length on the phase of linac LI is found to be in 
good agreement with the simulation predictions. Two cases of measured and simulated bunch time- 
profile are presented in Fig. ||. Again, the agreement between simulation and experiment is fairly 
good taking into account the uncertainties associated to the zero-phasing method. Noteworthy is 
the achieved peak current of ~150 A. 

The frequency-domain technique is based on the measurement of the coherent radiation emitted by 
the electron bunch via some electromagnetic process. In the coherent regime (i.e. for frequencies 
co ~ 2-7r/<7j where at is the rms bunch duration) the radiated power scales with the squared charge 
and depends on the bunch form factor. Thus it provides indirect informations on the bunch time- 
profile. In DUV-FEL, we detect the far-field radiation associated to the geometric wake field caused 
by aperture variation along the accelerator (e.g. the irises of the rf-structure). The radiation shining 
out of a single-crystal quartz vacuum window, located prior to the linac section L3, was detected 
with a He-cooled bolometer. The detection system, composed of the bolometer and the vacuum 
extraction port, can transmit radiation within the frequency range [cj/,u; w ] ~ [1.2,50] THz. The 
lower and upper frequency limits being respectively due to diffraction effects related to the finite 
size of the detector and transmission function of the vacuum extraction port. Given the bunch 
charge Q and the Fourier transform of the bunch time-profile S(oj), the power is expected to scale 



as P ~ Q 2 f£ u cLu\S(lu)\ 2 oc Q 2 /at (see annex for details). The typical signal observed as the charge 
is varied is presented in Fig. [?]: the observed nonlinear behavior confirms that the emitted radiation 
is not incoherent. From simulation we expect the power to scale as P oc Q 1,37 (see annex for details) 
a number close to the one resulting from the fit of the data: P oc Q 1 ' 57 . 

In Figure [8|, the measured bolometer output signal versus the phase of LI is compared with the 
expectation (1) calculated from the simulated phase space density and (2) computed from the mea- 
sured bunch time profile obtained by zero-phasing. As expected the increase of the coherent signal 
is an unambiguous signature of the bunch being compressed (the charge was monitored during the 
measurement and remained constant to 200±20 pC). 

The data points computed from the measured time profile fmeas were obtained by numerically com- 
puting the Fourier transform of the bunch time-profile (using a FFT algorithm) and by performing 
the integration: 

fmeas = / duj\S(uj) X R{oj)\ 2 . (8) 

where R(u) stands for the frequency response of the detection system. 

To generate the data points from the simulated phase space distributions f s imu we write 
the time-profile, S(t) as a Klimontovitch distribution: 

1 N 

i=l 

iV being the number of macro-particle used (50000 in the simulations presented hereafter) and t, 
the time of arrival of the i-th. macro-particle. Eq. (pf) allows to write the integrated power as: 
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(10) 



Though Figure § shows the signal increases as the bunch is compressed, there are discrepancies 
between the measurement and the two calculations for the short bunch case, we believe this is due 
to the lack of a precise knowledge of the transmission line frequency response. 



IV. CONCLUSION 



We have measured the bunch length dependence on the phase of a traveling wave accelerating 
structure located just downstream from an rf electron source. We could compress the bunch by a 
factor ^3, down to ~0.5 ps, for a bunch charge of 200 pC. In our experimental setup, a stronger 
compression is currently difficult to achieve without significantly impinging the transverse phases- 
pace quality. The linac section used for the compression also plays a crucial role in achieving low 
emittance since it quickly accelerates the beam at energies of approximately 60 MeV thereby freezing 
the transverse phase space. Hence operating the first linac far off-crest reduces the final energy 
and impact the emittance since transverse space charge forces scale as I/7 2 . An improvement of 
our experiment would be to surround the accelerating structure used as a bunch compressor with 
a solenoidal lens to enable a better control of the beam transverse envelope and emittance |0, 15|. 
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Appendix: Dependence of radiated power on bunch charge 



Let's consider the case of a Gaussian distribution: 



The corresponding bunch form factor takes the form: 

/+oo . . 

S(t) exp — iu)t\ 2 = exp (-a 2 uj 2 ) , 
-co 

and the integrated bunch form factor in the [w;,u; u ] frequency interval is: 

f = dco exp ( —a t lo ) . 

Ju>i 

The integration of the latter equation can be written in term of "error" function: 



/ = 77— (erf(crtW M ) - erf(a t uJi)) . (11) 



Taking into account the limit of the erf function, 



lim erf(z) = 1, and (12) 

z — >oo 

lim erf(z) = 0, 
z—>0 



and assuming the frequency range is so that otoj u 3> 1 and atuJi < 1, we finally have for the radiated 
power: 

O 2 

PocQ 2 x/oc— . (13) 
at 

Figure ^ shows the dependence of the bunch length versus the charge expected from simulations . 
We find at oc Q 0A3 and thus we would expect the radiated power to be P oc Q 1 ' 57 which is close to 
the value deduced from the fit of the measurement presented in Fig. [7|: P oc Q L37 . 
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FIG. 1: Simulation, using a two macro-particles model, of the velocity compression in a 3 m long traveling 
wave structure. The initial conditions are (i/> , 7o)=(0,8) and the macro-particle spacing is (Ai/> , A7 o )=(0.1, 
0.4). Plot a) shows snapshots at different z of the longitudinal phase space each segment extremities 
is determined by the two macro-particles positions. Plots b) and c) present the energy gain and phase 
evolution of the two macro-particles versus z. In these two latter plots, solid lines represent the leading 
particle and dashed lines the trailing one. 



klystron 1 



rf-gim 



klystron 2 klystron 3 







LI 




L2 





L3 




L4 







to 

undulator 



bunch compressor 



UV laser 



YaG monitor 



FIG. 2: Overview of the Deep ultra-violet free-electron laser (DUV-FEL) accelerator. LI, L2, L3, and L4 
are the four linac sections. 




FIG. 3: Simulation of the zero-phasing method for a distorted incoming longitudinal phase space. The 
images a) and c) depict the phase spaces after the bunch as passed the zero-phasing traveling wave structure 
for the "positive" (upper plots) and "negative" (lower plots) zero-crossing points. The plots b) and d) are 
the corresponding projections. In these plots we compare the time projection (dashed lines) with the one 
deduced from the fractional momentum spread projection (solid lines). The time > corresponds to the 
bunch tail. 




FIG. 4: Total energy versus phase of linac section LI. The points are measurements obtained parasitically 
to the bunch length measurement. The solid line is a simulation result. 




FIG. 5: Rms bunch length versus phase of the linac section LI. 





FIG. 6: Comparison of the bunch time-profile for LI on crest {4>li = 0°), and —82° off-crest. Plot a) was 
generated by tracking simulation; plot b) is a direct measurement using the zero-phasing method. The time 
> corresponds to the bunch tail 
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FIG. 7: Dependence of bolometer signal versus bunch charge. The circles are measurement, the line is a fit 
of the measurement using & a x Q@ law, the result gives (3 = 1.37 ± 0.06. 
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FIG. 8: Integrated bunch form factor / normalized to its value at <f)Li = 38°. fboiom, fmeas, and fsimu 
correspond respectively to measurement with the bolometer, computation from the measured time-profiles 
and computation from the simulation-generated time profiles. 




FIG. 9: Simulated dependence of bunch length versus charge. The circles are simulation results, the line is 
a fit to the simulation using a, a x law, the result gives (3 = 0.437 ± 0.007. 



